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ABSTRACT 

The one particle sector of the simplest one dimensional quantum lattice gas automaton has 
been observed to simulate both the (relativistic) Dirac and (nonrelativistic) Schrodinger 
equations, in different continuum limits. By analyzing the discrete analogues of plane 
waves in this sector we find conserved quantities corresponding to energy and momentum. 
We show that the Klein paradox obtains so that in some regimes the model must be 
considered to be relativistic and the negative energy modes interpreted as positive energy 
modes of antiparticles. With a formally similar approach — the Bethe ansatz — we find the 
evolution eigenfunctions in the two particle sector of the quantum lattice gas automaton 
and conclude by discussing consequences of these calculations and their extension to more 
particles, additional velocities, and higher dimensions. 
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1. Introduction 

The recent interest in quantum lattice gas automata (QLGA) is at least partly stimulated 
by the prospect of quantum computation. Quantum computers are (as yet hypothetical) 
physical systems designed and programmed to take advantage of the quantum mechanical 
phenomena of superposition and interference in order to perform calculations more effi- 
ciently than is possible with a classical computer [1] . As the size of semiconductor devices 
decreases, useful quantum coherence first seems possible at the nanoscale [2] where the 
most plausible computer architecture is that of a cellular automaton [3] . Thus it is natural 
to consider quantum computation with lattice gas automata (LGA). 

Just as classical LGA are more efficiently deployed to simulate physical processes like 
diffusion and fluid flow than for universal computation [4], it seems likely that, as antici- 
pated already by Feynman [5] , the first useful quantum computation will be a QLGA simu- 
lation of some quantum mechanical process. Thus it is important to understand the physics 
of QLGA: their symmetries, their conserved quantities, and their macroscopic/ continuum 
limits. 

In the next section we begin by reviewing the one particle sector of the simplest QLGA 
in one dimension. Determining the discrete analogue of plane waves, we find a complete set 
of invariants. These conserved quantities bear interpretation as energy and momentum, 
and enable the exact computation of the evolution. 

This model has been observed to simulate both the (nonrelativistic) Schrodinger equa- 
tion [6,7] and the (relativistic) Dirac equation [8,9,10]. Although perhaps surprising at first, 
this is not inconsistent as the Dirac equation is well known to have a nonrelativistic limit 
[11] and, as Benzi and Succi have emphasized, numerical evolution of parabolic equations 
can be stabilized and localized by the use of artificially hyperbolic equations [6]. In Section 
3 we use the physical interpretation of the conserved quantities in the model to distinguish 
relativistic and nonrelativistic regimes of the QLGA. Specifically, we show that the Klein 
paradox obtains [12]: a sufficiently discontinuous geometry /potential precludes interpre- 
tation of the QLGA as simulating a single nonrelativistic particle, compelling recognition 
of antiparticles in the model as well. 

In Section 4 we extend the one particle model to a multiparticle QLGA. Introduction 
of all the two particle and antiparticle scattering amplitudes makes the QLGA formally 
(almost) equivalent to the six vertex model [13]. The exact solvability of the six vertex 
model leads us to expect a computable complete set of invariants for the multiparticle 
QLGA. Concentrating on the two particle sector of the model, we find these invariants 
explicitly and explain their physical meaning. 

In Section 5 we conclude with a brief discussion of our results and their extension to 
more particles, additional velocities, and higher dimensions. 



2 



QLGA and their invariants 



David A. Meyer 



2. The one particle sector 

A LGA should be envisioned as a collection of particles moving synchronously from vertex 
to vertex on a fixed graph (lattice) L: At the beginning of each timestep each particle is 
located at some vertex and is labelled with a 'velocity' indicating along which edge incident 
to that vertex it will move during the 'advection' half of the timestep. After moving along 
the designated edge to the next vertex, in the 'scattering' half of the timestep the particles 
at each vertex interact according to some rule which assigns new 'velocity' labels to each. 
Here we will consider only the one dimensional lattice isomorphic to the integer lattice Z. 

A QLGA is a LGA for which the time evolution is unitary. To make this precise 
we must first identify the Hilbert space of the theory. An orthonormal basis for the one 
particle sector Hi of the Hilbert space H of the one dimensional QLGA is given by {\x, a)} 
(in the standard Dirac notation [14]), where x G L denotes position and a G {±1} denotes 
'velocity'. At each time the state of the QLGA is described by a state vector in Hi. 

x,ct 

where the amplitudes i/j a (t, x) G C and the norm of ^(t), as measured by the inner product 
on H, is: 

l=^p^ji(t,4 (2.1) 

x,ct 

The state vector evolves unitarily, i.e., ^{t + 1) = U^>(t), where U is a unitary operator 
on H. Since the evolution is unitary, the inner product is preserved and (2.1) holds for all 
times if it holds for one; this allows the interpretation of ip a (t, x)ip a (t, x) as the probability 
that the particle be in the state \x, a) at time t [14,15]. The total probability is the most 
fundamental conserved quantity of a QLGA. 



The 'advection' followed by 'scatter- 
ing' evolution of the quantum particle is 
implemented by having the basis vectors 
evolve as 

U\x, a) = a\x+a, a)+b\x+a, — a). (2.2) 




Figure 1. Spacetime diagrams of the one parti- 
cle scattering amplitudes. Time runs upward. 



That is, an initially right (left) moving quantum particle moves one lattice point to the 
right (left) and then continues in the same direction with amplitude a or changes direction 
with amplitude b as shown in Figure 1. We showed in [10] that the most general unitary 
evolution for a one dimensional QLGA with parity invariance* and particle speed 1 is 
described in the one particle sector by a unitary scattering matrix 



I.e., invariance under x — > —x; also called reflection invariance. 
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S may be parameterized by a = cos 9, b = ism 9. The model is then equivalent to 
Feynman's path integral for a Dirac particle of mass 9 [8], but one might ask for the origin 
of the evolution rule (2.2) and (2.3) — from what does the particle scatter? One answer to 
this question will emerge at the end of the next section. 



It is convenient to introduce the par- 
ity transformed scattering matrix S := 
PS, where 



1 

1 



b a 
— > 



a b 
<— 



x-1 



x+1 



Then (2.2) becomes 



Figure 2. Schematic diagram of a subspace of 
H±. Each box represents one dimension of the 
Hilbert space. The arrows indicate the amplitude 
flow into the two components at x; the labels a 
and b are the multiplicative factors to the corre- 
sponding boxes. 



U\x, a) = S a ' a \x + a, a') 



(2.4) 



This is equivalent to a condition on the amplitudes: 

Uip a/ (x) = ^2S a/a ifj a (x - a), 



(2.5) 



where the time argument is omitted. Figure 2 illustrates the amplitude flow in Hi described 
by (2.5). The discrete analogues of plane waves are wave functions if) a : L — > C which 
evolve by phase multiplication. Since U is unitary, these are exactly its eigenvectors; they 
evolve by multiplication by the corresponding norm 1 eigenvalues e~ luJ . For plane waves, 
therefore, (2.5) becomes 



e M Ipg' (x) = ^ Sgiqlpa (x-a). 



(2.6) 



Solving an equation like (2.6) is a problem we will face three times in this paper. In each 
case we proceed the same way, by making an appropriate ansatz for the form of a solution 
and then constraining the parameters in the ansatz by imposing additional (matching) 
conditions. 

As we explain in some detail in [16] for a QLGA which generalizes the one defined by 
(2.2) and (2.3), the appropriate ansatz for (2.6) is 



^ a k \x-a) = e- iak ^ a k \x), 
just as we would expect for a plane wave. Plugging (2.7) into (2.6) gives 



(2.7) 



(2.8) 
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As the notation suggests, the ip a {&) are the components of a vector ip(x) G C 2 (in this 
notation the state vector is ^ 
problem: 



i/j(x)\x)); thus (2.8) can be rewritten as an eigenvalue 
■-^^ k \x) = D(k)^ k) (x), 



where 



D(k) :-- 



ae lk 
be ik 



be lk 
ae~ ik 



is unitary since S is. D(k) summarizes the amplitude flow into the central pair of blocks 
in Figure 2. 



The characteristic equation for D(k) 
is the dispersion relation: 

cos uj = cos 9 cos k. (2-9) 

For pairs (u, k) satisfying the dispersion 
relation, the plane wave ansatz satisfies 
(2.6). The dispersion relation (2.9) has 
two solutions ±u;fc for each k, as shown 
in Figure 3. The eigenvalues e^ %UJk have 
eigenvectors 

-be~ ik 



CO 



^,(*.±)(0) 



ae 



(2.10) 



The plane wave solutions to (2.6) are thus 

^ k >±\x) = e lfcx V (fc,±) (0), (2.11) 



71/2 




-71 


-7C/2 




71/2 


K 






-7C/2 










-7C 







Figure 3. The dispersion relation for 9 = n/12. 
There are two branches and a frequency gap be- 
tween 9 and — 9. 



where for positive cu/c, positive (negative) 
k gives a right (left) moving plane wave. 
That the wave functions (2.11) are the dis- 
crete analogues of continuum plane waves is visible in Figures 4 and 5 which show, respec- 
tively, the evolution of right moving k = n/16 and k = n/8 plane waves for < x < 31. 




Figure 4. Evolution of the k = n/16 right mov- 
ing plane wave for the interval < x < 31 when 
9 = 71-/12. 



Figure 5. The same simulation for the k = ir/8 
right moving plane wave. This wave has a higher 
frequency. 
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Just as in the continuum situation, therefore, uj and k bear interpretation as being 
proportional to energy and momentum, respectively [16]. Furthermore, they label a com- 
plete set of conserved quantities in terms of which the evolution can be exactly computed: 
Let 

\k,e) := J2^ {k ' £) (°) elkX \ x ) 

X 

denote the plane wave state vectors in H\. Then (dealing with normalization in the usual 
way [8]) the set {|/c,e)} is an orthonormal basis for Hi [16]. Consider any state vector 
# ( //;: 

X 

= ^2^(x) \k,e)(k,e\x) 

x k,e 
k,e x 

The parenthesized expression gives the components in the new \k,e) basis: 

$e{k) := ^2(k,e\x)ip(x) 

X 

= ^(X,^ e) (0)e iky \y))\x)^(x) 

x y 

= (^ (fc ' e) (0)) t ^V(^)e" lfcx 

X 

=: (V (fc ' e) (0))V(*), 

where ip(k) is the discrete Fourier transform of ip(x). The plane waves \k,e) evolve by 

phase multiplication so the probabilities ifj e (k)ifj e (k) are left invariant by the evolution. 
Since any initial state vector \I/ (0) can be expressed in the plane wave basis this way, the 
existence of these conserved quantities is equivalent to exact solvability for the one particle 
sector of the QLGA. 

There is, however, a familiar problem with the physical interpretation of these con- 
served quantities. We have identified uj as being proportional to energy, so for each k there 
are plane waves with either positive or negative energy. We could restrict our state vectors 
to have nonzero amplitudes only in the directions \k, +), i.e., not allow any negative energy 
states; the conservation laws would preserve this condition. This would, however, preclude 
highly localized states: \1/ = \xq, +1), to take an extreme example, has nonzero amplitudes 
in all the negative energy directions \k, — ). Furthermore, introducing any interaction into 
the model risks disturbing this symmetry and allowing the system access to the (appar- 
ently) unphysical negative energy states. In the next section we will see exactly how this 
can happen and how we should reinterpret the model as a consequence. 
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3. The Klein paradox 

A single electron moving ballistically in a nanoscale semiconducting device is a conceiv- 
able implementation of a QLGA. Alternatively, simulation of such a physical situation is 
a possible application of a QLGA implemented in some other way. In either case, inhomo- 
geneous potentials must be included in the model. As explained in [17,8,7,16], modifying 
the evolution rule by an inhomogeneous phase factor has the desired effect. For a time 
independent potential, (2.5) is modified as 

Ui> a .{x) = Y J So t > a e- l * {x - a H a {x-a), (3.1) 

a 

where <j>(x) is the potential. 

Consider a simple potential step: 

and let us study right moving plane waves scattering off the step. Using (3.1), the equation 
analogous to (2.6) is 

e-^Mx) = J2^ a e- i(t,{x - a) Mx - a). (3.3) 

a 

For x < 0, inserting the potential (3.2) into (3.3) gives exactly equation (2.6). With 
an incident right moving plane wave we expect a reflected left moving plane wave, so an 
appropriate ansatz is 

^){x) = tfj (k ' +) (x) + A^~ k ' + \x) if x < 0, (3.4) 

where (a;, ±k) satisfy the dispersion relation (2.9). 

For x > 1, inserting the potential (3.2) into (3.3) gives 

e- i "il; a ,{x) = Y,S a >*e- i( t ) il> a {x-a) 

a 

which is equivalent to 

e-^-^ix) = Y,Sa'aMx-a). (3.5) 

a 

Again (3.5) has the same form as (2.6) so we make the ansatz 

^\x) =B^ k '^+\x) if x>l, (3.6) 
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b a 



-1 



a cb 







Figure 6. Schematic diagram of a subspace of H\ when there is a potential step of height 
<fi between x = and x = 1. As in Figure 2, the boxes represent dimensions of the Hilbert 
space and the labelled arrows indicate the amplitude flow (here c := e~ t ^). Equations (3.8) 
correspond to the statements that the left moving amplitude flow from x = 1 equals what it 
would be were the ansatz for x < continued to x = 1, and that the right moving amplitude 
flow from x = equals what it would be were the ansatz for x > 1 continued to x = 0, 
respectively. 



for the transmitted amplitudes, where 

cos(c<; — </>) = cos cos k'. (3.7) 



As illustrated in Figure 6, at x = and a; = 1 (3.3) gives a pair of boundary conditions 
which constrain the two parameters A and B: 

e --^)( ) = ^,_ ie -^y (1) + 
e--V# } (l) = S^-ie-^y (2) + ^+iV#?(0). 
These constraints imply that 

5e-^V-i'' +) (l) = V>_i +) (1) + v4Vi"i fc ' +) (l) 
^+) (0 ) + A^- fc ' +) (0) = S e -^V?i' +) (0), 



(3.8) 



after inserting the ansatz (3.4) and (3.6). Using the formulae (2.10) and (2.11) for the 
plane waves and simplifying, the pair of equations (3.8) becomes 

A-Be~ l4> = -1 
Ae~ lk - Be~ l ^- k,) = -e ik . 

These equations can be solved for A and B to give 

A= ' ' 



^/ik^ g %k 



B = e^ € ' 



^%k^ g ik 
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Figure 7. The real part of the a = — 1 com- 
ponent of an eigenfunction of U with frequency 
uj = 7r/6 for a step of height 7r/24 with = 7r/12 
on the interval — 128 < a; < 128. This wave 
function has wave numbers aproximately 0.459 
to the left of, and 0.296 to the right of the step. 




Figure 8. The eigenfunction of U in the same 
situation as in Figure 7 except that in this case 
the step is taller — it has height 7r/8. The incident 
right moving wave still has wave number approx- 
imately 0.459. To the left of the step the wave 
function decays exponentially. 



This justifies the ansatz of (3.4) and (3.6) as a solution to (3.3). 

Notice that the solution k' to (3.7) will not always be real, so (3.6) and (3.7) will 
not always describe a plane wave. For example, given < 9 < uj, when the step is small, 
i.e., < (f) < uj — 9, the transmitted amplitudes form a plane wave of lower frequency 
uj — (f) (and smaller wave number k') than the incident plane wave, but still lying on the 
upper branch of the dispersion relation shown in Figure 3. When the step is larger, i.e., 
uj — 9<4><uj + 9,uj — (p lies in the gap between the branches of the dispersion relation 
so (3.7) has purely imaginary solutions k' and the only physical solution of the form (3.6) 
decays exponentially past the step. 

Figures 7 and 8 illustrate these situations on the interval —128 < x < 128 with 9 = 
7r/12 for an incident right moving plane wave with uj = tt/6 and k = arccos(cosu;/ cos 9) ~ 
0.459. In Figure 7 the step has height 7r/24 and we have plotted the real part of V^i • 
This is the case Q<(f) = Tx/2A<uj — # = 7r/12 so there is a transmitted plane wave with 
k = arccos(cos(u; — (f))/ cos 9) ~ 0.296. In Figure 8 the step has height 7r/8; this is the case 
uj — 9 = n/12 < <j) = n/8 < uj + 9 = n/4 so there is no transmitted plane wave and the 
amplitude decays exponentially in x. 

This much is exactly what we would expect in the nonrelativistic situation; and we 
would expect the transmitted amplitudes to simply decay exponentially faster as the height 
of the step is increased. Instead, we find that when the step is even larger, i.e., uj+9 < 0, the 
transmitted amplitudes again form a plane wave, but now with frequency uj — <p < —9 < 0. 
Figure 9 illustrates this situation — the Klein paradox [12] — with a square well of depth 
(p = 7n/24 and uj = n/6 again, so that uj + 9 = n/4<(p = 7n/24. Here the transmitted 
plane wave has frequency uj — <fi = —tv/8 < and hence negative energy. This situation is 
out of the nonrelativistic regime; the large potential step localizes the particle too strongly, 
exciting negative energy states. 
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Just as in the continuum model for 
a relativistic quantum particle, therefore, 
since we take negative energies to be un- 
physical, we are forced to recognize the ex- 
istence of ant/particles in the model, i.e., 
particles which are the time reverse T (or 
charge conjugate C) of the particles with 
which we started. The tokens represent- 
ing the antiparticles must be the 'empty' 
(lattice point, direction) pairs; that is, the 
single particle state \x,a) should now be 
interpreted as a particle at (x, a) and an- 
tiparticles at all (y,(3) ^ (x, a). Since the 
antiparticles are the time reverse of parti- 
cles, a negative frequency corresponds to 
a positive energy state for them; thus there are no unphysical negative energy states. In 
this relativistic interpretation, the amplitude for the single particle to change direction is 
identified with the amplitude for it to scatter back from an oppositely moving antiparticle— 
answering the question asked in the previous section. A complete reinterpretation of the 
model requires definition of the two antiparticle, and hence two particle, scattering ampli- 
tudes. We do so, completing the definition of the QLGA, in the next section. 

4. The two particle sector 

In Section 2 we defined only the one particle sector of the model; in [10] we showed that 
the most general parity invariant local evolution rule which preserves particle number and 
has an exclusion principle is defined by the unitary scattering matrix 

( d \ 

a b v ' 

\ f J 

where the central block is the scattering matrix S (2.3) and unitarity requires dd = 1 = //. 
The phase / is the amplitude for a pair of particles to scatter off one another, while the 
phase d is the amplitude for either a pair of empty states to evolve to a pair of empty states 
or a pair of antiparticles to evolve to a pair of antiparticles, depending on our interpretation 
of the model. Figure 10 shows all the scattering rules implied by R; as in Figure 1 time 
runs upward so we indicate the presence of antiparticles by downward pointing arrows. 
This notation displays the (formal) equivalence of our one dimensional QLGA with two 
copies of the six vertex model [13]. The fact that two copies of the six vertex model are 
required stems from the usual Z2 symmetry in one dimensional LGA: particles an even 
distance apart initially cannot interact on the integer lattice [18,19]. 




Figure 9. The eigenfunction of U in the same 
situation as in Figures 7 and 8, but now for an 
even taller step, of height 77r/24. Rather than 
decaying even faster than in Figure 8, the ampli- 
tude of the wave function oscillates to the right of 
the step, with wave number approximately 0.296. 



R 
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* X * * X 

Figure 10. Spacetime diagrams of all the QLGA scattering amplitudes. Time runs upward so 
antiparticles are represented by downward pointing arrows. 

In the nonrelativistic interpretation it is most natural to set d = 1, while in the 
relativistic interpretation we must set d = f to maintain CT invariance [20,19]. In either 
case, the evolution rule (2.2) for the basis vectors of the one particle sector of the Hilbert 
space gets multiplied by a factor of d for each antiparticle/hole. Since an overall phase is 
unobservable, we can always multiply R by d to set Rn = 1 and avoid the normalization 
problem associated with an infinite number of factors of d. The one particle sector of 
the model is now completely defined and so supports the particle/antiparticle relativistic 
interpretation introduced at the end of the last section. 

Furthermore, since particle number is preserved, only the relative phase between S and 
/ has observable consequences, so the parameterization of S in the previous sections can 
be maintained in each of the multiparticle sectors of the QLGA. Since the inhomogeneous 
potentials considered in Section 3 are also implemented by a relative phase, this form for the 
scattering matrix R suggests that similar calculations will determine the eigenfunctions of 
the evolution matrix for interacting particles. In this section we demonstrate this explicitly 
for the two particle sector of the QLGA. 

An orthonormal basis for the two particle sector H2 C H is given by the set { \x\, cti)® 
\x2,ol-2) I {x\iOL\) 7^ (^2 5 «2)}; this reflects an exclusion principle in the model — the two 
particles cannot occupy the same state — but not necessarily fermionic statistics. A state 
vector for a pair of particles has the form 

ai«2 

(t,xi,X2)\x 1 ,ai) <g> \x2,a 2 ) 

Oi,ai)/(x 2 ,a 2 ) 

and just as in the one particle sector must have norm 1: 



ai«2 

(t,X!,X2). 

The evolution of the two particle basis vectors implied by the scattering matrix R 
(4.1) divides into two cases. For x\ + ol\ 7^ X2 + 0.2 the evolution rule (2.4) is simply 
applied to the two particles independently: 

U^ax) <g) \x2,a 2 ) = £ S a > iai S a > 2a2 \x 1 +a 1 ,a' 1 )®\x2 + a 2 ,a2), (4.2) 
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aa ba 
ab bb 



(x+1, x-1) 



Figure 11. Schematic diagram of a subspace of if2,int- The coordinates x\ and x-i of the two 
particles label the sets of four boxes, each of which corresponds to a pair of 'velocities'. The 
evolution rule (4.4) gives the amplitude flows shown pointing toward the diagonal (where the 
particles coincide), while (4.5) gives the amplitude flows shown pointing away from it. In order 
for the ansatz (4.9) to satisfy (4.7) for x\ ^ X2, the amplitude flows away from the diagonal 
must satisfy (4.10); in order for it to satisfy (4.11) the amplitude flows toward the diagonal 
must satisfy (4.12). 



while for x\ + a\ = x 2 + cx 2 = x, there are two possible scattering events: 

U\x - 1, +1) <g> \ x + 1, -1) = f\x, +1) <g> \x, -1) 
U\x + 1, -1) ® \x - 1, +1) = f\x, -1) ® \x, +1) 
which can be summarized as 

U\x — a, a) <S> \x + a, —a) = f\x, a) <E> |x, —a). (4-3) 

The evolution rules (4.2) and (4.3) can be rewritten in terms of the amplitudes: For 

xi ^ x 2 , 

U^ a ' lCe ' 2 (x 1 ,X 2 ) = ^ Sa'^Sa^a^a^ix! - a U X 2 ~ a 2 ) (4.4) 
ai ,a 2 

and for x\ = x 2 = x, 

Uif) a - a (x,x) = fip a - a (x -a,x + a). (4.5) 

Just as Figures 2 and 6 illustrate the evolution rules (2.5) and (3.1), respectively, Figure 
11 illustrates the amplitude flow on H 2 described by (4.4) and (4.5). 
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Notice that we have implicitly assumed that the two particles are distinguishable, as 
would be the case classically, and have interpreted the scattering events with amplitude 
/ to involve neither particle changing direction. One might be concerned that (4.3) and 
(4.5) would have to be changed for indistinguishable quantum particle scattering. This is 
not the case when the particles have fermionic statistics, for example: if the initial wave 
function is antisymmetric under the exchange of two particles, i.e., 

1p aia2 (x U X 2 ) = -1pa 2 a 1 (x2,Xl), (4.6) 

then the evolution rules (4.4) and (4.5) preserve this property. We will consider this case 
in more detail at the end of this section. 

At the beginning of this section we noted the close relation of our one dimensional 
QLGA with the six vertex model. The six vertex model is exactly solvable [21,13], which 
suggests that we should be able to find the analogue of plane waves in the two particle 
sector, namely the two particle eigenfunctions of U. The associated conserved quantities 
will describe the physics of this sector. For an eigenfunction, (4.4) becomes 

e~ luJ ip a ' 1 c t ' 2 {xi,x 2 ) = ^2 ^c.!^^^!^^! -<*i,x 2 -a 2 ). (4.7) 

ai ,a 2 

A straightforward ansatz for a solution to (4.7) is: 

^ ai a 2 (x llX2 ) = V£ 1 ' £l) (^l)V'£ 2 ' e2) (^), (4.8) 

where -i/A fci ' £i ) =: ip^ are the one particle plane waves we found in Section 2. It is easy to 
verify that (4.8) solves (4.7): 

(x u x 2 )= ^^^a/" 1 " 1 " 1 ^^!)^^^!^) 

= E^' 1 « 1 e "^ lfcl < ) (^)E^^ 2 e " ia2fc2 < ) (^) 

= e~ l ^iji 1 }(x 1 )e-^^(x 2 ) 
= e -^+^^ a , ol ,(x 1 ,x 2 ), 

where the third equality follows from (2.8) since the pairs (e^, ki) satisfy the dispersion 
relation (2.9). Thus the ansatz (4.8) solves (4.7) for uj = eiUi + e 2 u 2 . 

Now observe that the Z 2 symmetry of the one dimensional QLGA partitions H 2 into 
two orthogonal subspaces: H 2 = H 2jint © #2, free, where 

#2,int := span{ \x 1 , a{) <g> \x 2 , a 2 ) \ x 1 - x 2 = (mod 2) and (xi, cti) 7^ (x 2 , a 2 ) } 
H 2 j Tee ■= span{ \xi, cti) <g> \x 2 , ol 2 ) \ X\ — x 2 = 1 (mod 2) }. 
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These subspaces evolve independently since H 2y i nt consists of states in which the two 
particles can eventually interact while #2, free consists of states in which they cannot, as 
is implied by the form of (4.2) and (4.3). In fact, the evolution of H 2 ^ ree is completely 
specified by (4.2) (or (4.4)) and depends only on the one particle scattering matrix S. On 
this subspace, therefore, (4.8) provides a complete solution for the eigenfunctions of U. 

To find a solution for the interacting particles subspace H 2 ^ nt , consider the two particle 
plane wave (4.8) 'incident' from the x\ < X2 subspace of i?2,int- This two particle plane 
wave 'scatters' off the x\ = xi subspace of i?2,int just as the right moving plane waves we 
considered in Section 3 scattered off the potential step, producing both a 'reflected' and a 
'transmitted' two particle plane wave, in the x\ < X2 and x\ > x 2 subspaces, respectively. 
Thus we make the Bethe ansatz [22,21,13]: 

,/,(w) „ \ _ J ^alM^cS M + AlpS (x 2 )lp ( ah x l) if x l < x 2] (A q\ 

[ Bf^(x 1 )f^{x 2 ) if Xi > X 2 , 

where u> := e\uj\ + e 2 u 2 , since all three two particle plane waves on the right hand side 
of (4.9) have the same two wave numbers h\ and k 2 and the same frequency 10. For 
l^i — CC2 1 7^ 2, the free particle subspace computation together with linearity imply that 
(4.7) is satisfied by (4.9). For \x\ — x 2 \ = 2, (4.7) implies conditions on the extension of 
this ansatz to x\ = x 2 = x: 



^€U X - !) = E *«i«Ai« a <«> -l-a u x + l-a 2 ) 

(a 1 ,a 2 )#(-l, + l) 

+ S a ^ ^S^^l^x.x) 
e- l ^l"U x + 1^ " 1) = E S a , a S a> /:^{x + l- ai ,x-l-a 2 ) 

(ai,a 2 )^(+l,-l) 

+ ^ 'a'^+l^ ^,-1^+1,-1 x ) ■ 

These conditions will be satisfied provided 

^-l+ifrx) = ^\{x)^\{x) + A^\{x)^\{x) 

(4.10) 

^i- 1 {x,x) = B^i{x)^i{x), 

i.e., provided the form of the ansatz (4.9) extends to the two nonzero components of 
^L 2 (x,i) in this way. 

Now the problem is to find A,BeC such that (4.9) and (4.10) define an eigenfunction 
of (4.5); that is, must also satisfy 

e-^^l^x, x) = f^ljx -a,x + a). (4.11) 
14 
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This gives two constraints on A and B: 

e"*" [^l(x)^l(x) + A^l(x)^l(x)} = f Br^lix + l)^\{x - 1) 

e-^B^l (a;) v£? (x) = / - (x + 1) + ^V'-i (z + Ity+i (ar - 1)] . ' " 



Inserting the form (2.11) of the one particle plane waves into (4.12) and rearranging gives 

(4.13) 

where ip eie2 := ipi? {Q)^ (0). Solving (4.13) for A and B gives 



Ae i(fc 1 -fc 2 )^_ + _ Be -™^ = _ e -<(fci-fc 2 )^ + _ 



A= V-+^+-(/ 2 -e- 2 ^ 



5 



(e-^+_) 2 - (e i ( fc i- fc 2)/ ? /;_ + )2 
(e"^ + _) 2 - (e l ( fc i- fc 2)/ ? /;_ + )2 



(4.14) 



Thus (4.9) and (4.10), with the coefficients A and S given by (4.14), define a set of 
eigenf unctions for the two interacting particles subspace H 2 ^ nt . Each is described by a pair 
of wave numbers k\ and k 2 , and a frequency u; = eicui + £2^2, where the pairs (eja>j,A;j) 
satisfy the dispersion relation (2.9). When the two particles are apart, i.e., not at the 
same lattice point, these eigenfunctions each put the particles into a superposition of two 
particle plane waves, with momenta proportional to k\ and k 2 . Thus we should interpret 
these two particle eigenstates as describing pairs of particles with distinct momenta before 
they interact (i.e., arrive at the same lattice point), and the same distinct, but possibly 
permuted, momenta after they interact. The relative amplitudes of the permutations are 
A and B in (4.14). To understand their physical meaning, notice that \A\ 2 + \B\ 2 = 1. 
Since A and B generically have both nonzero norms and phase angles, one may interpret 
B, for example, as describing the relative probability of transmission |£?| 2 as well as a 
phase angle shift 0(ki, k 2 ) between the incoming and scattered (transmitted) two particle 
plane waves. 

When A ^ 0, a necessarily independent set of eigenfunctions for the two interacting 
distinguishable particles subspace H 2jini is obtained by the same calculation, only starting 
with an 'incident' two particle plane wave in the x\ > X2 subspace: 

iptl (xi,x 2 ) = { B '^K x i)^}{x 2 ) if (xi,ai) < (x 2 ,a 2 ); ( 415 ) 

where the ordering on pairs (x^, ctj) is lexicographic and A' and B' are given by the formulae 
(4.14) for A and B, but with k\ and k 2 exchanged (including the k\ and k 2 in ip |_ and 

V>+_). 
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These eigenfunctions (4.15), together with our first set (4.9), form a basis in terms of 
which antisymmetric wave functions (4.6) (and symmetric ones) can be expanded. Since 
indistinguishable fermionic particles with a necessarily antisymmetric wave function are 
the most plausible components of a QLGA quantum computer, however, it is perhaps 
more natural to work in the antisymmetric subspace of H 2 from the beginning. Then the 
natural eigenfunction ansatz is antisymmetric: 



V'iV \xi)Wa} \x 2 ) + A^a}{x 2 )^a}{xi) if (xi,tti) < (x 2 ,Oi 2 ); 
-^a}{xi)^a}{x 2 ) - A^{x 2 )^a}{xi) if fa, Oil) > (X 2 , Oi 2 ). 

(4.16) 

The ansatz (4.16) satisfies (4.7); the interaction condition (4.11) imposes only a single 
constraint in this case: 

Inserting the form (2.11) of the one particle plane waves and solving for A gives 

A e-^-+ + fe- i ^- k ^+- 

Thus (4.16) with the coefficient A given by (4.17) defines a set of antisymmetric eigenfunc- 
tions in the two interacting particles subspace; the eigenfunctions in the two free particles 
subspace are simply obtained by antisymmetrizing (4.8). 



5. Discussion 



We began our analysis in Section 2 by determining the eigenvalues and eigenfunctions — the 
phases e~ tuJ and the plane waves ^( fc ' e ) = \k,e) — of the evolution operator U on the one 
particle sector of the QLGA. We found that the probability that a single particle with 
initial wave function ^(0) has wave number ko and frequency ujq = eo arccos(cos 9 cos ko) 
at time t is an invariant: | (^(0) \ko, eo) | 2 - It is constant in time since each plane wave \k, e) 
evolves by phase multiplication e~ lujt over t timesteps. Consequently, the expectation value 
of the wave number 

(k) :=$>|<vl>(0)|fc,e>| 2 

k,e 

and of the frequency 

(u;}:=J2^\^(0)\k,e}\ 2 

k,e 

are also invariants of the evolution. Physically, of course, this means that momentum and 
energy are conserved. 

Identification of k and u> as proportional to momentum and energy, respectively, led 
to our interpretation in Section 3 of the eigenfunctions of U in the presence of a large 
step potential not as negative energy modes of the particle but rather as positive energy 
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modes of antiparticles in the model. Since there are no antiparticles in the (nonrelativistic) 
Schrodinger equation continuum limit of this model, our analysis of invariants in the one 
particle sector of the QLGA enables us to identify this as a parameter regime whose 
continuum limit must be described by the (relativistic) Dirac equation. 

The eigenfunction ansatz and matching conditions we used for the step potential 
forshadowed the Bethe ansatz approach we applied in Section 4 to find the eigenvalues 
and eigenfunctions of U in the two particle sector of the QLGA. Just as the one particle 
plane waves determine evolution invariants in the one particle sector, the two particle 
eigenfunctions \k\, e\\ kz, 62) '■= tp^ defined by (4.9) and (4.10), or (4.15), and (4.14) (or 
by (4.16) and (4.17)) define invariants in the two particle sector: The probability that a 
pair of particles with initial wave function \I/ (0) have wave numbers ki and frequencies 
Ui = Ci arccos(cos 9 cos ki) at time t is an invariant: |(\l/(0)|/ci, ei; /?2, £2)| 2 - Again, this is 
simply because the two particle eigenfunction evolves by phase multiplication by e~ tuJ , 
where to = e\U)\ + €2^2- 

In fact, the Bethe ansatz extends to each of the multiparticle sectors of the QLGA 
[22,21,13] so the evolution of an arbitrary initial wave function can be computed exactly. 
Of course, it had to be the case that the unitary evolution operator U has eigenfunctions 
which evolve by phase multiplication; what is remarkable, and what is meant by 'exactly 
solvable', is that these eigenfunctions and eigenvalues can be effectively computed, just as 
we have done in the one and two particle sectors of the QLGA. 

The most general, no less local, QLGA in one dimension includes particles with 've- 
locity' [10]. We have determined the one particle plane waves for this model in [16] and 
it is natural to ask if the Bethe ansatz also solves this model exactly. One expects this 
calculation to be possible, but more difficult: H2, for example, does not decompose into a 
direct sum of 'interacting' and 'free' subspaces as it does in the single speed case analyzed 
here. Similarly, Boghosian and Taylor have simulated the evolution of one particle plane 
waves in a two dimensional QLGA [7] . One would like to extend the Bethe ansatz to find 
multiparticle eigenfunctions in such a higher dimensional model, but should expect it to 
be even more difficult — there are very few higher dimensional lattice models known to be 
exactly solvable [13,23]. 
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